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Abstract 

We investigate the steady-state solution and it's bifurcations in time-delay systems 
with band-limited feedback. This is a first step in a rigorous study concerning the 
effects of AC-coupled components in nonlinear devices with time-delayed feedback. 
We show that the steady state is globally stable for small feedback gain and that 
local stability is lost, generically, through a Hopf bifurcation for larger feedback 
gain. We provide simple criteria that determine whether the Hopf bifurcation is 
supercritical or subcritical based on the knowledge of the first three terms in the 
Taylor-expansion of the nonlinearity. Furthermore, the presence of double-Hopf bi- 
furcations of the steady state is shown, which indicates possible quasiperiodic and 
chaotic dynamics in these systems. As a result of this investigation, we find that 
AC-coupling introduces fundamental differences to systems of Ikeda-type [Ikeda et 
al., Physica D 29 (1987) 223-235] already at the level of steady-state bifurcations, 
e.g. bifurcations exist in which limit cycles are created with periods other than the 
fundamental "period-2" mode found in Ikeda-type systems. 
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1 Introduction 



Recently there has been considerable interest in chaotic devices operating in 
the radio-frequency (RF) regime for applications such as ranging and commu- 
nication. As an example, such devices are studied as a possible signal source 
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for radar applications because chaotic waveforms have the desirable proper- 
ties of a large frequency-bandwidth and a fast decay of correlations [1,2]. Fur- 
thermore, microwave [3,4], optoelectronic [5,6,7,8], and optic [9] devices are 
considered for communication systems since they can generate chaos with fre- 
quencies that match the frequency range of the communication infrastructure 
and provide advantages such as increased privacy [7] and high power efficiency 
[4]- 

Many RF-devices are most accurately modeled by delay differential equations 
(DDEs) because the time it takes for signals to propagate through the device 
components is comparable to the time scale of the dynamics. Furthermore, 
many of the chaotic devices are designed explicitly to include a time-delayed 
feedback [3,4,5,6,7,8,9] because it introduces advantageous features such as the 
ability to tune the complexity of the dynamics by adjusting the delay [10] and 
the possibility of control in the presence of substantial control loop latency 
[11]. In modeling, it is important to take the frequency characteristics of the 
feedback into account. 

At high speed, many components are AC-coupled, which means that low- 
frequency signals below the frequency cut-off are suppressed. As a conse- 
quence, the time-delayed feedback signal is band-pass filtered because in ad- 
dition to the cut-off at low frequencies, high frequencies are suppressed due to 
the finite response time of device components. Thus, DDEs describing band- 
limited feedback are the appropriate model for many of the devices that have 
been developed. 

In the literature, we find cases where a model for the device-dynamics is con- 
structed which ignores high-pass filtering in the feedback, but the experimen- 
tal device has AC-coupled components. See for example, Ref. [3]. This suggest 
that the effect of AC-coupling on the dynamics is considered to be ignorable. 
However, counter-examples exist, such as the work of Goedgebuer et al. [12] 
and Blakely et al. [8]. Goedgebuer et al. find that the inclusion of a high-pass 
filter in the feedback drastically increases the dimensionality of the chaotic 
attractor when compared to the same system without such filter [12]. Simi- 
larly, in our work on a high-speed chaotic time-delay system [8] , we found that 
it was necessary to model the effects of AC-coupling in order to explain the 
observed bifurcation of the steady state and the route to chaos. These counter- 
examples demonstrate that high-pass filtering due to AC-coupled components 
can influence the dynamics substantially, indicating that a better theoreti- 
cal understanding of the effects of high-pass filters in time-delay systems is 
needed. As a first step in this direction, this paper provides a detailed mathe- 
matical analysis of the bifurcations from the steady-state solution of systems 
with band-limited feedback. 

The essential building blocks of the class of delay systems studied in this 
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Fig. 1. Sketch of a time-delay system with passive nonlinearity G, amplifier of gain 
7, and band-limited delayed feedback. 

paper include a passive nonlinear element, a feedback loop, and a way to 
provide amplification (see Fig. 1). The nonlinear element maps the input to 
the output, x out = G(xi n ), the feedback loop is responsible for the delay T, 
and the control parameter 7 is a measure of the signal amplification. Taking 
into account the band-pass characteristics of the feedback results in models 
that are two-dimensional DDEs in the simplest case, where, for simplicity, it 
is assumed the transfer characteristics of the feedback can be approximated 
by a two-pole band-pass filter. 



dzi(t) = zi(t) + dz 2 ^ 
dt r h dt ' 



n- 



dz 2 (t) 
dt 



-z 2 {t)+ 1 G Zl {t-T) 



Here, the time-scale is related to the corner frequency of the high-pass filter 
through uj h = t^ 1 and r\ is related to the corner frequency of the low-pass filter 
through uji = 7] _1 . 

In this paper we study both the global and local stability of the steady-state 
solution of (1) for general nonlinear functions G. We derive explicit equations 
for the Hopf-bifurcation curves and use center-manifold techniques to pro- 
vide simple criteria for the type of these bifurcations. Furthermore, we show 
that double-Hopf bifurcations of the steady state are possible, which indicates 
possible quasiperiodic and chaotic dynamics in these systems. 

To the best of our knowledge, Hopf bifurcations of (1) have not been studied in 
detail. On the other hand, the scalar DDE that result by taking the limit where 
the corner frequency of the high-pass filter goes to zero has been studied in- 
tensely both experimentally [13,14,15] and theoretically [16,17,18,19,20,21,22], 
starting in 1979 with the pioneering work of Ikeda [16]. In the discussion sec- 
tion we will therefore compare delay systems with band-limited feedback and 
Ikeda-type scalar DDEs to pinpoint the new dynamic features that arise due 
to the presence of the high-pass filter. 

The paper is organized as follows. In Sec. 2 we present the model equations. 
In Sec. 3, a sufficient condition for global stability of the steady-state solution 
is given. In Sec. 4, equations determining the boundary of local stability are 
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derived and it is shown that the steady state becomes generically unstable 
via a Hopf bifurcation. The main focus of this section is the solution of the 
relevant characteristic equation. In Sec. 5, center-manifold theory is used to 
determine whether the Hopf bifurcation is subcritical or supercritical. In Sec. 6, 
three numerical examples are discussed. We conclude with a discussion section, 
Sec. 7. 



2 The Model 



To simplify our analysis, we bring model (1) in a more convenient form by 
introducing variables t, x, and y through 

t = t (- + -), x = z u y ^^^( Zl - Z2 + 1 G[0}), (2) 
\n T h J r h + Ti 



to obtain 



<h ' [f) --x(t)+y(t) +1 f[x(t-r)], 



dt 

dy(t) 
dt 



= —rx(t). 



(3) 



Here, the dimensionless delay is r = T (r t 1 + r h 1 ), r = TiT h {ji + r h ) 2 , and 
f[ x (t - r)] = r h (r h + Ti)-\G[x(t - r)] - G[0]). The nonlinear function / is 
defined such that /(0) = 0. 

There are three dimensionless parameters that influence the dynamics: the 
gain 7, the strictly positive delay r, and r (0 < r < 1/4). The parameter 
r is related to the angular frequency at which the transfer-function of the 
bandpass-filter is maximum. Indeed, the frequency that maximizes transmis- 
sion is u m ax = \Juiu)h (see Fig. 1), which, in the new coordinates, corresponds 
to the dimensionless angular frequency 

i tmax = ; = Vr. (4) 

Additionally, r provides a measure of the bandpass-filter bandwidth, because 
one may uniquely identify small r with a large bandwidth and r < 1/4 with 
a narrow bandwidth if one assumes that Uh < Ui holds. 

The trivial steady-state solution, i.e. x = y = 0, is the only fixed point of 
model (3). The local stability of the fixed point can be studied for general 
nonlinear functions by considering the Taylor-expansion of f(x) around x — 0. 
Retaining the lowest-order term of the Taylor-expansion and multiplying by 
the gain factor 7 yields the effective slope, one of the relevant bifurcation 
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parameters. Let us therefore introduce 



b = 7 /'(0). 



(5) 



The value of the effective slope b where the fixed point first becomes linearly 
unstable is called the critical value, denoted as be- It determines the critical 
gain through ^ c — f'i^)~ l bc- 

For the discussion of linear stability and for the center-manifold analysis it 
is useful to rewrite (3) as a functional differential equation and to explicitly 
separate the linear and nonlinear parts. Model (3) may be written as 



u{t) = L ll {u t ) + F{u t ,n) (i>0), 



(6) 



where we use the following notation. Let u(t) G R 2 be the vector u = (x,y) T , 
and let C — C ([— r, 0], R 2 ) be the Banach space of continuous functions map- 
ping the interval [— r, 0] (r > 0) into R 2 . For (p G C the norm is defined as 
\\ip\\ = sup_ T<e<0 | |, where | • | is the norm in R 2 . Furthermore, for each 
fixed t, u t designates the function in C given by u t {6) = u(t + 6), 6 G [— r, 0]. 
The parameter ji is defined through fj, = b — bc, i.e. it is the bifurcation param- 
eter shifted such that /i — at the critical value. L M : C — > R 2 is a bounded 
linear operator and F is a C k function containing only nonlinear terms, i.e. 
F(0,/i) = and the first derivative D 1 F(0,fi) = 0. The linear operator L M is 
given by 

LM = j_ r dr](0, t*M0) := L (<p) + ^(cp), (7) 

where 77 is a n x n matrix- valued function of bounded variation on [— r, 0], 
(p E C ([— r, 0], R 2 ), and 



, -111 b c 0\ 10' 

%)= U(0)+ U(-r), L^tp) = 

-r \00/ \00, 



<P(-r). (8) 



The function F is given by 







3! 



r) + /i.o.t. , (9) 



where y^ 1 -* denotes the first component of the vector p = (ip( 1 \ip^) T and 
the coefficients 62 (fJ>) and 63 (/i) are related to the Taylor-expansion of the 
nonlinearity /(#) through the relations 



W=7/ ff (0) = (6c + |i)^, 
W=7r(0) = (6o + |i)^. 



(10) 
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We will use the abstract description of time-delay systems with band-limited 
feedback given by Eq. (6)-(10) to determine the stability of the fixed point 
u = and to derive the Hopf normal-form. However, we will refer to the more 
concrete description provided by model (3) when we state results. This will 
make it easier to apply our findings to experiments. 



3 Global Stability of the fixed point for small gain 

In the discussion of stability properties of solutions to dynamical systems the 
notions of global and local stability have to be distinguished. Global stability 
implies that the fixed point is the only attractor of the system, whereas local 
stability means that, after small perturbations away from the fixed point, the 
system will return to it's original state. In this section, a sufficient condition for 
global stability of the steady-state solution of delay systems with band-limited 
feedback is provided. 

Consider the case 7 = (no feedback) reducing (3) to the second-order equa- 
tion of a damped harmonic oscillator. In this case, the fixed point is globally 
stable. For small amounts of feedback one would expect that the dissipation 
in the system will dominate and the fixed point will as well be the global 
attractor. Under some assumptions on / we can show that this is true (see 
Appendix A for the proof). 

Proposition 1 Consider (3) with f(x) being continuous, Lipschitz in x, and 
in particular assume that there exists a constant kf such that kf\x\ > \f{x)\ 
for all i: / 0. // |7| < kj 1 , the trivial solution of (3) is globally uniformly 
asymptotically stable. 

The assumption that / is Lipschitz is not very restrictive because most phys- 
ically implementable nonlinearities are differentiable (have a bounded deriva- 
tive) and are therefore Lipschitz. For these systems it is always possible to find 
a small but finite gain I7I for which the condition \x\ > |7||/(a;)| (Vx 7^ 0) of 
Proposition 1 is satisfied. However, in many cases the steady-state solution will 
be globally stable for a larger range of I7I because Proposition 1 only provides 
a sufficient condition. Furthermore, the bound was derived so that it holds for 
arbitrary delays r and therefore may not be optimal. Nevertheless, the bound 
is both sufficient and necessary for nonlinearities / for which kf is minimal, 
i.e. kf = |/'(0)|i and thus global stability extends to |6| = |7/'(0)| < 1. The 
reason Proposition 1 provides as well the necessary condition for these non- 
linearities is that, for \b\ > 1, there exist values of r for which the steady state 
is locally unstable as will be shown in the next section. 
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4 Local Stability of the Fixed Point 



To gain a better understanding of how different dynamics arise as system 
parameters are varied, the local stability of the fixed point is discussed in this 
section using linear stability analysis. The main idea is to investigate how small 
perturbations to the trivial solution evolve for a given set of parameter values, 
which is equivalent to studying the corresponding characteristic equation. The 
goal is to construct a bifurcation plot showing the location of bifurcations as 
parameters are varied. When a system has two or more parameters, as is 
the case here, it is common to draw two-parameter bifurcation plots. In this 
context, co dimension-one bifurcations are located on one-dimensional curves 
and codimension-two bifurcations occur at isolated points in the bifurcation 
plot. 

4-1 Characteristic Equation 

For (6), the local stability of the fixed point is determined by the eigenval- 
ues (or characteristic values) of the linearized equation, ii(t) = L^iut). The 
eigenvalues A satisfy the characteristic equation [23,24] 

det A(A, 0) = 0, A(A,//) = AI - L„ (e A l) . (11) 

and coincide with those values of A for which there is a nonzero vector c A such 
that there is a solution of u(t) = L^Uf) given by u t = c\exp(\t). With 
defined by (7), the characteristic equation (11) is given by 

A 2 + A + r - b\e- XT = 0. (12) 

To determine the local stability of the fixed point, we first consider solutions 
of (12), which is transcendental and has an infinite number of roots for every 
set of fixed parameter values. The fixed point is locally stable if all roots 
(eigenvalues) have a negative real part. Thus, to determine parameter values 
for which the fixed point becomes unstable, we set Re(A) = and Im(A) = iQ. 
Separating real and imaginary part yields 

= -fi 2 + r-6fisin(fir) 

= l-6cos(fir). ^ 

Note that this set of equations is unchanged for f2 — > — f2 and that there 
is no solution for Q = if r > 0. This implies that a pair (or pairs) of 
complex-conjugate eigenvalues will cross the imaginary axis as the bifurcation 
parameter is varied. Thus, it is sufficient to consider positive Q only and, 
generically, the trivial solution becomes unstable through a Hopf bifurcation, 
which is a co dimension-one bifurcation. 



7 



One way to visualize the solution of (13) is to seek parameterized curves in 
the plane of two bifurcation parameters, which we choose as the delay r and 
the effective slope b. These curves separate regions in parameter space with 
different numbers of eigenvalues in the right complex halfplane. The relevant 
region where the fixed point is stable (no eigenvalues in the right complex 
halfplane) is the one that includes 7 = 6 = 0. The parameterization of the 
curves is most conveniently achieved through s = Qt yielding 



rS(s) = — (tan(s) + ^Un'(s) + 4rj , (14) 

= Jw- (15) 

Here, the label n denotes the different solution branches and the value of n is 
defined in terms of the curve parametrization as 

0<s<! & n = 0, (16) 

(2n- 1)tt (2n+ 1)tt , n n 

^<s<^ & n = l,2,3,.... (17) 

Note that Vq is positive for even n and negative for odd n. Also, note that Tq 
approaches zero as s — > 0. Furthermore, Tq — * as s approaches (2n — l)7r/2 
from above but Tq tends to infinity as s approaches (2n + l)7r/2 from below. 

In Fig. 2, the stability boundary of the fixed point is shown in r-5-parameter 
space. The stability boundary is obtained by combining those pieces of the 
parameterized curves (tq(s), 6^(s)), given by (14) and (15), for which, upon 
crossing, the number of eigenvalues in the right complex halfplane changes 
from zero to two. From Fig. 2a it is seen that the fixed point is always locally 
stable if — 1 < b < 1, that it may be stable for a considerably larger range of 
b if r is small, and that \bc\ ~ 1 if t is large. 

Figure 2b provides a zoomed view of part of the parameter space where the 
parameterized curves (tq(s), b^{s)) are shown as solid lines (starting from the 
left the curve-index n is n = 2, 4, 6, 8, . . .). We will refer to these curves as 
Hopf-curves. As one of the Hopf-curves in Fig. 2b is crossed from below, a 
pair of complex conjugate eigenvalues moves into the right complex halfplane, 
which is seen by considering how the eigenvalues A change with respect to b 
for fixed r. Using (12) we calculate the real part of the derivative of A at the 
critical point, 



b n T C S (f + COS ( g ) Sin ( g )) + 25 2 COS 2 (s) 



critical 



\2s + rgtan(s) + t^scos(s) l e 



Because s = Qt > implies that s + cos(s) sin(s) > 0, it follows that the right 
hand side of (18) is always nonzero and that it's sign is determined by the 
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Fig. 2. Local stability of the fixed point, indicated by the shaded region, as a function 
of the parameters r and b. a) A large part of the stability region (— 1 < b < 1) has 
been contracted to make details of the boundary visible, b) The inset shows a portion 
of the stability boundary. The number of eigenvalues in the right complex halfplane 
is given. The square symbols mark points on the stability boundary where two 
pairs of complex conjugate eigenvalues lie on the imaginary axis, whereas crosses 
denote the extrema of the Hopf-curves, i.e. 6 = 1. The dashed lines depict the 
approximation given by (20). 

sign of Vq. Therefore, for branches with positive b^, the eigenvalues move into 
the right complex halfplane as b is increased, whereas they move into the right 
complex halfplane for branches with negative Vq as b becomes more negative. 

The above discussion confirms that the steady state will, generically, become 
unstable through a Hopf bifurcation in systems with band-limited feedback. 
Let us make this statement more precise. The two conditions for a Hopf bi- 
furcation are that there is a pair of simple characteristic roots crossing the 
imaginary axis transversally at /i — (b — be) and that the characteristic 
equation A(A, 0) has no other roots with zero real parts. The roots (eigenval- 
ues) cross transversally because the right hand side of (18) is always nonzero. 
Furthermore, there is exactly one pair of roots on the imaginary axis for all 
parameter combinations of b and r that fall on one of the Hopf-curves. The 
exception are points where two Hopf-curves intersect, because in that case two 
pairs of eigenvalues cross into the right halfplane and a codimension-two bi- 
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furcation (double- Hopf bifurcation) occurs. Let us denote as r + (t~) the set of 
delays r at which the fixed point becomes unstable due to such a codimension- 
two bifurcation as the effective slope b is increased (decreased). Because the 
sets r + and r~ are countably infinite, i.e. have Lebesgue measure zero on M + , 
the fixed point loses stability through a Hopf-bifurcation for essentially all 
delays r. In other words, careful tuning of the gain (oc b) and the delay (t) 
is necessary to bring an experiment close to a codimension-two bifurcation 
point. It is in this sense that the Hopf bifurcation is generic for time-delay 
systems with band-limited feedback. 

4-2 Approximate analytic solution of the characteristic equation 

Although (14) and (15) allow us to determine numerically the critical param- 
eters to a desired precision, it is useful to have an explicit expressions for the 
critical value of the bifurcation parameter. For many experiments the feedback 
gain is most easily varied, which means that the transcendental equations (14) 
and (15) should be inverted to give fe c as a function of the remaining two pa- 
rameters r and r. 

An exact solution is not possible in general. Therefore, we solve approximately 
around the extrema of the stability boundary, i.e. around points along the 
Hopf-curves where s = nrr (crosses in Fig. 2): 

777T" 

T £(s = mr) = — , b n c {s = mr) = ±1, fi c (s = mr) = yfi. (19) 

We expand the trigonometric functions and the square root in (14) and (15) 
around the extrema and solve for the critical effective slope 

/ 2r 2 (V-^) 2 2(wrr 5 / 2 + 4r 3 ) (V - 

b n c (r,r) = ± 1 + — -A %- - ; V V~r) +h o f (2Q) 

CK ! \ (n7r + 2v^) 2 (n7r + 2v^) 4 I v ; 

The approximation given by (20) is shown in Fig. 2 using dashed lines. 
4-3 The frequency at the onset of instability 

So far we have discussed the values of the feedback gain for which the steady 
state is stable. That is, we have shown how to calculate the critical gain j c f° r 
known values of the delay (r), the first derivative of the nonlinearity (/'(0)), 
and the frequency of maximal transmission (yjr). The value of 7c can be 
compared with experiments. Another useful way to compare experiments and 
theory is to study the frequency of oscillations at the onset of instability. 
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Fig. 3. Imaginary part of the eigenvalue versus the delay, ^(t) is shown for positive 
b (n is even) and r = 0.028 . The part of each branch n that belongs to the stability 
boundary is shown as bold straight line. The dashed lines are the approximate value 
of the frequency. 



Consider the case where the Hopf bifurcation is supercritical. Then, as the 
stability boundary is crossed, the fixed point becomes unstable, a stable limit 
cycle is born, and the system starts to oscillate. Experimentally, the frequency 
of this oscillation is the quantity that is most readily measured. The frequency 
at onset is determined by Vt^ through = n^(27r)~ 1 (^ + ooi). Therefore, 
it is useful to plot the eigenvalue V£q versus the delay r as is done in Fig. 3. 



In the vicinity of the extrema (Eq. (19), crosses in Fig. 2), the frequency scales 
roughly as ~ nir/r (f£ ~ n/(2T)) for n > 0. This approximate behavior, 
shown as dashed curves in Figure 3, can be understood by considering whether 
a wave circulating in the feedback loop will reinforce itself. For the case of 
positive feedback (b > 0) a periodic perturbation will reinforce itself, if the 
feedback delay is a multiple of the wave's period, i.e. / ~ n/(2T) with n an 
even integer. On the other hand, a sinusoidal perturbation is amplified by 
negative feedback (b < 0), if it is shifted by half it's period after one round- 
trip. Thus, for b c < the frequency is expected to scale as / ~ n/(2T) with 
n an odd integer. This reasoning is consistent with the fact that for n even 
(odd) the critical effective slope b^ is positive (negative). 

From Fig. 3, it is also seen that different oscillation "modes" will be observed 
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as t is increased. That is, the frequency of the observed oscillations will jump 
from Q ~ nir/r to Q ~ (n + 2)h/t as r is varied across one of the double- 
Hopf points (r G t^; squares in Fig. 2 inset). These jumps are explained by 
the fact that the gain in the feedback loop is not perfectly flat over the pass- 
band. As 7 is increased from a low level, one particular frequency will first 
reach the threshold where the gain in the loop balances the losses. In a system 
with only low-pass feedback, the gain is highest at low frequencies, so the 
oscillation-mode with the lowest frequency is always the one that destabilizes 
the steady state, independent of the delay. On the other hand, the high-pass 
filter introduces a bias toward high frequencies. Because the frequency scales 
roughly as Q ~ r _1 for each mode n, the damping effect of the high pass filter 
on a particular mode becomes more pronounced with increasing delay time r. 
Therefore, there exists a delay r for which a higher order mode, one that has a 
higher frequency for a given delay, will reach threshold first. Thus, examining 
the frequency of the oscillations at the onset of instability as a function of the 
delay allows one to distinguish time-delay systems with band-limited feedback 
from systems with low-pass feedback [8]. 

For completeness, note that fl^ can be determined exactly if r is considered 
as bifurcation parameter and b is kept constant, and is given by 

(nn? = f I {[(b) 2 + 2r - 1] + v^) 2 + 2r - 1]* - 4r^} ifr£<^, 
U {[( fe ) 2 + 2r - 1] - v/[(&) 2 + 2r - 1]* - 4r*} if r£ > f r , 

as is discussed in detail in [25]. 



5 Hopf Normal-Form 



In the previous sections, we have shown that a time-delay system with band- 
limited feedback undergoes Hopf bifurcations as the system parameters are 
varied. However, we have not yet determined whether the Hopf bifurcations are 
subcritical or supercritical. The expected dynamical behavior close to thresh- 
old is very distinct for the two bifurcation types. For a supercritical Hopf bi- 
furcation, one expects small amplitude sinusoidal oscillations past the critical 
gain, whereas for a subcritical bifurcation, one expects bistability and hystere- 
sis close to the bifurcation point. This section is devoted to the determination 
of the type of the Hopf bifurcation using normal form theory. 

Normal forms have been studied extensively for finite-dimensional ordinary 
differential equations (ODEs) [26]. Recently, Faria and Magalhaes developed 
normal-form theory for the case of DDEs [27,28]. As for ODEs, one tries to 
simplify the dynamical systems under study using two approaches: reduce the 
dimension of the system and remove as much of the nonlinearity as possible. 
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The first goal is achieved by considering the restriction of the flow to the 
center manifold. This restriction is justified because the qualitative dynamical 
behavior of solutions of the infinite-dimensional DDE close to the fixed point 
can be described by the finite-dimensional ODE associated with the flow on 
the center manifold. The second goal is achieved by transforming the nonlinear 
differential equation into an equation with a simpler algebraic form, called a 
normal form, by changes of variables that eliminate all irrelevant terms from 
the equation, but retain the qualitative properties of the flow. 

The normal form theory developed in Refs. [27,28] uses the formal adjoint 
theory for linear DDEs [23] to set up an appropriate coordinate system near 
an equilibrium point. It then proceeds by performing a sequence of transfor- 
mations of variables such that, at each step j, the change of variables effects 
simultaneously the projection of the original DDE on the center manifold and 
removes the non- resonant terms of order j from the ODE on it. 

For a Hopf bifurcation, the normal form on the center manifold is given in 
polar coordinates p, £ by 

p = p K lP + K 2 p 3 + 0(p 2 p + | (jm, p) | 4 ) 
i = -S% + 0(\(j*,p)\). 

Here, p = b — is the bifurcation parameter, fl^ is the imaginary part of 
the eigenvalue, and K\ and K 2 are real coefficients. The qualitative behavior 
of the asymptotic solutions of (6) (or equivalently of (3)) is the same as the 
behavior of solutions of (22), which, in turn, only depends on the signs of the 
two coefficients K x and K 2 . Of special interest is to distinguish between the 
cases K 2 < indicating a supercritical bifurcation, and K 2 > indicating a 
subcritical bifurcation. The main object of this section is to express Ki and 
K 2 in terms of the relevant parameters (b, r, r, and the Taylor-expansion 
coefficients of the nonlinearity F) and to determine the sign of K\ and K 2 . 

5. 1 Derivation of the normal form coefficients 

For scalar DDEs (delay systems with low-pass feedback) an explicit formula 
for K\ and K 2 is given in Ref. [27] and reproduced in [24] (8.3). Based on this 
work, Giannakopoulos and Zapp [20] discuss in detail the different possible 
scenarios of Hopf bifurcations in scalar DDEs. As far as we are aware, no 
analogous work exists for two-dimensional DDEs. Nevertheless, we find that 
the derivation of the formula for K\ and K 2 for the two-dimensional case 
carries through entirely analogous to the derivation given in Refs. [24,27]. We 
will therefore restrict ourselves to the definition of the necessary quantities and 
to the statement of the end result. We adopt a notation that follows closely 
that of Ref. [24]. 
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We take //, or equivalently the effective slope b, as our bifurcation parameter. 
The other parameters, e.g. r and r, are assumed to be fixed. For a given set 
of parameters r and r, there are an infinite number of critical effective slopes 
distinguished here by the solution-branch label n (see (16) and (17)). For 
each be, there exists a center manifold at u — and the dynamics of (6), 
restricted to this manifold, is described by (22). That is, the derivation of the 
normal form is the same for all n. Furthermore, all expansion coefficients are 
evaluated at the bifurcation point, /i — 0. Consequently, to simplify notation, 
we will suppress explicit reference to n and ji for the rest of the section and 
write b := r := r£, fi := fig,, 6 2 := ^(a* — 0), and b 3 := b 3 (fx = 0). 

The first step in applying center-manifold theory to our problem is to introduce 
suitable local coordinates. For the Hopf-bifurcation it is convenient to work 
with complex variables, i.e to work in C = C ([— r, 0],C 2 ). Considering the 
linearized flow 

ii{t) = L (u t ) (t > 0), (23) 

one can associate with each fixed set of eigenvalues A = {Ai, . . . , \ s } a gener- 
alized eigenspace P with base <3> consisting of the appropriate eigenfunctions. 
In our case A = {iQ, — ifi}, and, denoting the complex conjugate using an 
overline, the base $ can be chosen as <3> = (0i, 2 ) with <f> 2 = 4>i, where 



i i 



ir/fl 



e iQe with 0! G C ([-r, 0], C 2 ) , -r < 6 < 0. (24) 



Using this notation, it follows that P = span$. To compute projections onto 
this generalized eigenspace one needs to consider the adjoint equation (see [23] 
for a definition). Let C 2 * be the space of row 2- vectors, C* := ([0, r], C 2 *), and 
let P* denote the dual of P. The base of P* is ^ = col(ipi,ip2) with ifj 2 = ipi 
and 

M0=n(i, -i/fi)e- inc with ^eC*([0,r],C* 2 ), < C < r. (25) 
The normalization factor is chosen to be 

N=(l + bre-' mr + V 1 . (26) 



fi 2 

The formal duality is the bilinear form (-|-) from C* x C to the scalar field 
[23,24] 

(alp) = a(0M0) - f f a{i - 6) drj(6) <p(t) a6C',^ C, (27) 
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(28) 



The normalization factor N (26) is chosen such that (^|$) = I, where I is 
the identity matrix. The purpose of above definitions is that it allows one to 
decompose the phase space C into the generalized eigenspace of (23) associated 
with A, i.e. P, and Q = {<p e C : {^\<p) = 0}. In short, the phase space C is 
decomposed by A as C = P © Q. 

In Ref . [28] , it is shown that the phase space associated with the full nonlinear 
problem (6) may similarly be decomposed and a projection from that phase 
space to P is given in terms of the bases $ and Furthermore, the center 
manifold theorem assures that there is a 2-dimensional invariant manifold of 
(6) tangent to the center space P of (23) at zero. It is therefore possible to 
use this projection and appropriate coordinate transformations to obtain an 
ODE on the center manifold in normal form, i.e. (22). We refer the reader to 
Ref. [24], Sec. 8.3 for details. 

Since the Hopf bifurcation is determined generically up to third order (see 
(22)) we need the Taylor-expansion of the nonlinearity F up to third order. 
Write the Taylor expansion of F at the critical point as 



Let H 2 : C x C — > C 2 be the bilinear symmetric form such that F(u, /i) = 



where the definition of the nonlinearity F (Eq. (9)) is used. 

Denote the complex coordinates on the two-dimensional center manifold by 
z, and the components of z by z\ and z 2 . These coordinates are related to the 
polar coordinates of the normal form (22) through z\ = pcos£ — ipsin£ and 
z 2 = zi. In the computation of the coefficients K\ and K 2 , the expansion of 
F(<J»z, 0) plays an important role. For the second-order terms write 



The coefficients are obtained by using the definition of (24), and F, (9): 




(29) 



F(u,0) = H 2 (u,u), 



(30) 



F 2 ($Z, 0) = 4(2,0,0)2? + 4(1,1,0)^1^2 + 4(0,2,0)^2- 



(31) 



4(2,0,0) = H 2 (<j)i,(j>i) = (o)b 2 e 




(32) 



4(1,1,0) = 2# 2 (0i,0 2 ) = (l)2b 2 



(33) 
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Write the third-order terms as 

F s (*z,im)= £ A (qi ^ l)Z fzf^ 1 . (34) 

<?l+<32+/=3 

For the Hopf bifurcation, the only relevant third-order coefficient is 

A( 2> i,o) = (5) 36 8 e- i0r . (35) 

With this result, we are finally able to state the main result of this section. 
Our calculation of K 1 yields 

K 1 = Re (^(0)^(00), (36) 

and for K 2 we obtain 

K 2 = I Re (^i(0)j4 (a>1>0 )) + \ Re {Vi(0)# 2 b&o (I)]" 1 Ai,i,o)) } 

+ ^ife{^l(0)if 2 (^, 

The above expressions are quite general [29] and are an obvious extension of 
the formulas for the Hopf normal-form coefficients for scalar DDEs given in 
[24], pp. 154, Eq. (8.77). 

Next, we express K\ and K 2 explicitly in terms of the parameters that are 
specific to the problem at hand. First, note that the following identity holds 

(38) 

critical 

Thus, the coefficient K\ is given by (18), which implies that K\ is positive 
for b > and negative for b < (see section Sec. 4). Second, in (37) the use 
of the definition of L , (8), and of H 2 , (30), and the substitution of ^4( 2 ,i,o), 
-4(2,0,0), 4(i,i,o), 0i, and Vi (given, respectively, by (32), (33), (35), (24), and 
(25)) yields 



2ifil - L (e 2im i) 



A (2,o,o)e 



(37) 



K x = Re 



~db 



Ko 



-Re| 
+ Im 



(brn 2 + (r + n 2 )e iQr ) / 



+ 



(b 2 ) 2 Q 3 



(brQ 2 + (r + tt 2 )e inT )(2ittb + {4Q 2 - 2iQ - r )e 2 ^) 



(39) 
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5.2 The sign of the normal-form coefficient K 2 



Equation (39) can be used directly to determine numerically the sign of K 2 
because the values of r, r, Q, and b on each of the Hopf curves are known. 
However, since we are only concerned with the sign of K 2 it is possible to 
derive a somewhat simpler criterion for distinguishing subcriticality from su- 
percriticality. To that avail, define a function C n (r, r) through 



C n (r,r) = y , (40) 

V 1 (b 2 tt 2 r + r + tt 2 ) \2ittb + (AVI 2 -2itt-r)e 2in T\ 2 

N n = [r + Q 2 ] [b 4 + 26 3 + b 2 - 4] + 2Q 2 rb 2 [b 3 - 1] 

+sgn(r - ^)3fiv / ^T[(4 - b 2 ) (r + Q 2 ) + 2Q 2 rb 2 ] . (41) 

Here, we restored the explicit reference to the label n for clarity. Recall that 
we use b = b^ and Q = Qq in this section and that the numerical values of 
b and Q depend on the remaining parameters r and r. It can be shown (see 
Appendix B) that 

Proposition 2 For r,r G M + ; r < \, and n — 0, 1, 2, . . . : 

• lff'(0) f"(0) + / /, (0) 2 C n (r,r) > 0, the Hopf bifurcation is subcritical (K 2 > 
0). 

• If /'(°) /'"(°) + /"(0) 2 C' n (r,r) < 0, toe #op/ bifurcation is supercritical 

(k 2 < o;. 

Since the Taylor-series coefficients of / are known, the difficulty in determining 
the sign of K 2 is shifted to finding the value of C n (r, r). Two simple cases that 
arise in this context are the following: 

Corollary 3 

(1) The case when there are no quadratic terms in the Taylor expansion of the 
nonlmearity (f"(0) = and /'"(0) ^ 0): If f (0)/'"(0) > 0, the Hopf bifur- 
cation is subcritical. If f'(0)f"'(0) < 0, the Hopf bifurcation is supercritical. 
In particular the type of bifurcation is independent of the sign of 7. 

(2) The case when there are no cubic terms in the Taylor expansion of the 
nonlinearity (f"'(0) = and /"(0) 7^ 0^: The Hopf bifurcation is subcritical 
ifC n {r,r) > and supercritical ifC n (r,r) < 0. 

In light of (2) in above Corollary, it is clearly useful to look for conditions 
on t and r for which the sign of C n (r, r) can be determined without explicit 
evaluation. To that avail we establish the following properties (see Appendix 
C): 
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Proposition 4 For r, r G M + ; r < \, and n — 0, 1, 2, . . . : 

• //r>^, C w (r,r)<0. 

• If r = ^ and n even {i.e. b = 1), C n (r, r) = 0. 
.If T< niLandbe (1,2), C n (r,r) > 0. 

• Ifr<^ and be [-2,-1], C n (r,r) < 0. 

Let us restrict our attention to the stability boundary of the fixed point, 
since this is most relevant for experiments. It is possible to show [30] that 
the stability boundary is bounded by 1 < b < 2 for all delays r (r > 0) for 
positive effective slopes. For negative effective slopes the stability boundary is 
bounded by — 2 < b < — 1 for r > t^(s) = r^(47r/3), where Tq(s) is given by 
(14). Combining this with the results of Proposition 2 and Proposition 4 the 
following simple cases arise in addition to those of Corollary 3: 

Corollary 5 At the stability boundary of the fixed point, the following holds 
for the Hopf bifurcation : 

(3) lfb<0, t > Tq{-^-), and /'(0)/'"(0) < ; the bifurcation is supercritical. 

(4) lfb>0, r > ^= 7 and f'(0)f"'(0) < 0, the bifurcation is supercritical. 

(5) lfb>0,r<^, and /'(0)/" / (0) > 0, the bifurcation is subcritical. 

In this section, we provide in Proposition 2 a criterion that determines the 
type of Hopf bifurcation in general but requires the numerical evaluation of the 
function C n (r, r). The Hopf bifurcation type can be determined without the 
evaluation of C n (r, r) if certain assumptions about the parameters of model (3) 
are satisfied. These simple cases are summarized in Corollary 3 and Corollary 
5. To illustrate the findings of this section we discuss in the next section some 
specific examples of time-delay feedback systems with band-limited feedback. 



6 Examples 

6.1 'Sine'-Nonlinearity 

As a first example, we consider f(x) = sin(x) as the nonlinearity in (3). In 
this case, there are no quadratic terms in the Taylor-expansion, f"(0) = 0. 
Therefore, according to Corollary 3, the Hopf bifurcation is always supercriti- 
cal because /'(0)/"'(0) < 0. This is confirmed by the numerical results shown 
in Fig. 4, which are obtained using DDE-BIFTOOL [31]. Displayed are the 
average amplitude of the periodic solutions versus the bifurcation parameter b, 
which is proportional to the feedback gain. Figure 4 also demonstrates that all 
bifurcations are supercritical not only for b > but as well for b = jf'(0) < 0, 
thereby confirming that supercriticality is independent of the sign of 7. 
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Effective Slope b Effective Slope b 



Fig. 4. The amplitude of the limit-cycle versus the effective slope b is shown for 
r = 2n/yfr and r = 0.028. a) As b is decreased past b = —1.009, the fixed point 
becomes unstable via a supercritical Hopf bifurcation (the Hopf curve with label 
n = 1 is crossed). As b is decreased further, additional limit-cycles are created via 
supercritical bifurcations as Hopf curves with label n = 3, 5, . . . are crossed, b) The 
fixed point is unstable for b > 1. All bifurcations are supercritical. 



Let us discuss in more detail the case of positive effective slopes b shown in 
Fig. 4b. As b is increased past 6=1 the fixed point becomes unstable and a 
stable limit cycle is created. A further increase of b results in the creation of 
additional limit cycles through supercritical bifurcations at the critical values 
of b. These additional limit cycles are stable within the center-manifold, but 
unstable with respect to the whole phase space, since the center-manifold itself 
has repelling normal directions. As these unstable limit cycles grow, additional 
bifurcations may occur that can change the stability property of the periodic 
solution. For example, in Fig. 4b, the periodic-solution branch emanating from 
b ~ 1.03 is unstable for small amplitudes, as predicted by our theory, but is 
stabilized through a torus-bifurcation (Neimark-Sacker bifurcation of the cor- 
responding map) at b ~ 1.07. In this context, it is important to emphasize that 
normal-form theory is valid sufficiently close to the fixed point and therefore 
makes statements only about small amplitude limit cycle solutions. 



Finally, the sufficient condition for global stability provided by Proposition 
1 yields the maximal range |6| < 1, since kf = f'(0) = 1, consistent with 
supercritical Hopf-bifurcations. 
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Effective Slope b Effective Slope b 

Fig. 5. a) Theoretical results: Shown are in the b — t— parameter plane the Hopf 
curves (6J(s), Tq(s)), given by (15) and (14) with r = 0.028. The curves correspond 
to the location of Hopf bifurcations. Dashed lines indicate subcritical and solid lines 
supercritical bifurcations, b) Numerical results: The periodic solutions that arise 
are explored for r = 27r/y / r ~ 37.5 (the dash dotted line in the left panel crossing 
Hopf curves with n = 2,4,6,8). Except for the n = 2 branch, which is a marginal 
case, the bifurcations are subcritical. 

6.2 Nonlinearity without Cubic Terms 



As a second example, we consider a nonlinearity without cubic terms in the 
Taylor-expansion 

f(x) = (x + -x 2 + x 3 )e~ x \ (42) 
2 

According to Corollary 3 the type of the Hopf bifurcation is determined by 
the sign of C n (r, r). Proposition 4 establishes this sign for all possible values 
of r and r under the condition that |6| < 2. 

To illustrate, we show in Fig. 5a the Hopf curves in the 6 — r— parameter plane 
and indicate whether the Hopf bifurcation on these curves is predicted to be 
subcritical or supercritical. 

Since the bifurcation is supercritical for r > nn / ' ^fr and subcritical for r < 
mr/y 7 ?, it follows that the case r = rm j 'yjr has to be marginal. Marginal 
means that the periodic solution branch emerges vertically when the amplitude 
of the limit cycles that are created at the bifurcation point is plotted versus 
the bifurcation parameter b. This is shown in Fig. 5b, where r ~ 2n / ' ^fr and 
it is seen that the periodic solution branch created from the n = 2 curve 
at b = 1 emerges essentially vertically. The other periodic solution branches 
shown in Fig. 5b are all found to be created via subcritical Hopf bifurcations, 
in agreement with the predictions of Corollary 3 and Proposition 4 (compare 
Fig. 5a and Fig. 5b). 
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6. 3 Experiment of Blakely et al. 



As our third and final example we consider the optoelectronic device described 
in Ref. [8]. We first map the model given in [8] onto (3) and then show that 
our theory correctly predicts the experimental results. 

In [8], the following model is derived for the optoelectronic device with band- 
pass feedback 

T h P(t) = - (P{t) - P ) + r h KV(t), 

n V(t) = - V(t) + 7 G [P(t-T D )] , (43) 



G 



P(t-T D )] =P(i-T D ){l + /3 sin [a (p(t - T D ) - P )]} , 



where P denotes the measured optical power and V the voltage in the elec- 
tronic feedback loop. The parameter values are taken from Table I in Ref. [8]. 
These equations can be cast in the form of (3) by introducing the rescaled 
and dimensionless variables t = t^r^ 1 + t^ 1 ), x = (P — Po)/Po, and y = 
T h (r h + n)- 1 {(P - P ) - k(V - 7G[0])} P \ The nonlinear delay term / of 
(3) is defined as f(x(t - r)) = r h (n + r^ 1 {G[P (x(t - r) + 1)] - G[P ]} P l 
and is given by 

f(x) = — (x + (3 x sin(ax) + [3 sin(ax)) . (44) 

Th + n 

The fixed parameters of the model are r = 29.8, r = 0.028, a = aP = 
49.14, r h = 22 ns, r { = 0.66 ns, and (5 = 0.8. For the Taylor-series of /, we 
obtain /'(0) = 39.14, f'(0) = 76.33, and f"'(0) = -92163. The values of the 
bifurcation parameter b attainable in the experiment are b G [0, 3.28], where b 
is defined through b := «ry/'(0). 

To illustrate the results of our analysis, we determine the critical gain and 
the type of the Hopf bifurcation using the above model. The boundary of the 
local stability region can be estimated using (20). We obtain be = 1.003, which 
corresponds to a gain 7 = 5.34 mV/mW, in agreement with the experimental 
result of 7 = 5.1 ± 0.5 mV/mW (see as well Fig. 3 of [8]). In the experiment, 
it was found that the bifurcation is supercritical. Since f'{0)f"'(0) < 0, it 
follows from Corollary 5 that the bifurcation is supercritical if b is negative 
(negative gain) and is supercritical for positive b if r > m\j\fr. However, a 
close inspection of the inset of Fig. 2 reveals that, for the given parameters (r = 
29.8 and r = 0.028), the steady state loses its stability as the parameterized 
curve with branch-label n = 2 is crossed to the left of the extrema, i.e. r = 
29.8 < 1-nj^fr. Consequently, their experiment does not fall into one of the 
simple cases, necessitating a numerical evaluation of C n (r, r), which reveals 
that /' (O)f'(O) + f"(0) 2 C n (r,r) < 0. Thus, the bifurcation is predicted to be 
supercritical in agreement with the experimental result. 
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7 Discussion 



In this paper, we take a first step in analyzing rigorously the dynamics of two- 
dimensional DDEs that model delay-systems with band-limited feedback. We 
provide a condition guaranteeing global stability of the steady state solution 
and use linear stability analysis to derive the boundary of the local stability 
region in parameter space. We show that, generically, the steady state will lose 
stability in a Hopf bifurcation. The use of center-manifold techniques allows 
us to give a criterion for the type of this bifurcation which holds for arbitrary 
nonlinear functions /. Some effort went into simplifying this criterion to make 
it easier to predict whether supercritical or subcritical behavior will occur in 
a given experiment. 

As mentioned in the introduction, the limit where high-pass filtering is re- 
moved from the feedback, i.e. in the limit r — > in (3), the two-dimensional 
DDE studied in this paper reduces to the class of scalar DDEs that includes 
the Ikeda DDE. It is informative to compare our results to known facts about 
Ikeda-type DDEs, since this brings into focus the novel features that arise due 
to the presence of the high-pass filter. 

The first and obvious difference is that blocking feedback at zero frequency 
in band-pass systems prevents the occurrence of multiple coexisting steady 
state solutions that exist for the case of scalar DDEs. Consistent with this 
observation, we find that fold-bifurcations do not exist in band-pass systems. 
For scalar as well as two-dimensional DDEs, Hopf bifurcations occur and give 
rise to limit cycles. However, there is a fundamental difference concerning the 
frequencies (or modes) of these oscillations. For Ikeda-type DDEs, the limit 
cycle at onset is always the fundamental "period-2" mode, with a frequency 
f ~ (4T)~ X for small delays and a frequency / ~ (2T) _1 for large delays 
(T ^> t{) [17,22]. In contrast, high-pass filtering results in a stability boundary 
where the mode at threshold varies with the chosen delay. As an example, in 
delay systems with negative band-limited feedback (b < 0), the n = 1 mode 
is observed for small delays. This mode corresponds to the "period-2" mode 
of Ikeda-type DDEs because it has a frequency / ~ (4T) _1 for T xs and 
a frequency / ~ (2T) _1 for T = 2nu;^ ax . As the delay is increased past the 
timescale associated with maximum transmission (T = 2noj^ ax ) ) modes with 
n — 3, 5, . . . (/ ~ n/(2T)) form the stability boundary. Thus, AC-coupling 
introduces fundamental differences to Ikeda-type DDEs already at the level of 
steady-state bifurcations, and the effects of high-pass filtering are particularly 
pronounced for large delays. 

Also in this paper, we show that double-Hopf bifurcation points exist for delay 
systems with band-limited feedback, whereas these cannot exist in first-order 
scalar DDEs with a single delay. However, since double-Hopf points are found 
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in n'^-order scalar DDEs [32], which may be used to represent multiple-pole 
low-pass-filter feedback, their appearance should be seen as a result of the 
increase in degrees of freedom rather than as a result of high-pass filtering. 

Let us conclude with some remarks concerning the double-Hopf bifurcation 
points. In experiments, the chance of choosing a fixed delay such that upon 
changing the gain the stability boundary will be crossed close to a double-Hopf 
point is negligible. However, the existence of double-Hopf points can provide 
insight concerning possible dynamics for feedback gains above threshold. As 
an example, double-Hopf interactions lead to limit cycles as well as tori [26]. 
Consistent with this observation, torus attractors were found numerically in 
[8] for the example in Sec. 6.3. Furthermore, chaotic dynamics exists near 
a generic double-Hopf bifurcation [26]. It would therefore be interesting to 
study how the coefficients of the double-Hopf normal form are related to the 
parameters of (3) and to determine whether there are restrictions on the flow. 
Restrictions means that the mapping of the parameters of the DDE onto the 
normal-form coefficients does not allow certain combinations of the double- 
Hopf normal-form coefficients to occur and thereby forbids the qualitative 
dynamics associated with those combinations. That such restrictions might 
occur was shown for instance by Buono et al. for the case of n 4/l -order scalar 
DDEs [32]. 



A Appendix: Proof of Proposition 1 

We are interested in investigating the stability of the trivial solution of the 
initial- value problem for (3), i.e. an equation of the form 

Z(t) = g(t,x t ), (A.l) 

where the operator g(t,(p),g : K. x B H — > R d is continuous and Lipschitz in 
tp G Bh and g(t,0) = 0. 

To prove Proposition 1, we use an extension of the method of Liapunov- 
Krasovskii functionals due to Kolmanovskii and Nosov [33]. With the Liapunov- 
Krasovskii approach, difficulties arise because it is necessary to construct a 
functional with negative-definite derivative. It is often easier to find Liapunov- 
functionals with nonpositive derivative. For this case stability may still be 
proved. 

To show stability of the trivial solution, two continuous functionals V(t, <p) and 
W(t, ip) are used. We require that the derivative of W is integrally unbounded. 
More precisely, designate Bh as the closed ball of radius H in the Banach space 
C = C([-r,0],R d ) 

B H = {<peC,d(0,<p)<H}, 
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where d(x,y) = \\x — y\\ denotes the distance. The the derivative W is called 
integrally unbounded in a set S C Bh if, for any number K > 0, there exists 
a number T(K) > and a continuous function £(t) such that uniformly in 
x t = x(t + 9) G S for t > t 

rt+T(K) 

We now cite the two theorems needed for the proof: 

Theorem 6 (Theorem 5.7, pp. 77, [33]) It is necessary and sufficient for 
uniform asymptotic stability of the trivial solution of (A.l) that there exist the 
functionals V(t, (p) and W(t, <p) such that 

(1) wi (|p(0)|) < V(t,<p) < u 2 (\\<p(6)\\), t > t , <p{ff) e B H , t e K. i/ere, 
CL?i, i = 1,2, are some scalar, continuous, nondecreasing functions, such that 
u>i(0) = and u>i(s) > for s > 0. 

(2) V" < £73(3^) < 0, where u73(y?) «s a continuous functional defined on Bh; 

(3) \W(t,<p)\ <L,t> t , ^(0) G BhM e 

(4) /or any /x G (0, if) there exists a p > snca taat tae derivative W is 
integrally unbounded in the set E{p, p) C C , where 

E(p, p) = M9) G B H , d (tp(9), S(ZJ 3 = 0)) < p, p < \\<p(e) \\ < H} . 

Here S(oJ 3 = 0) = & G B H ,oJ 3 ( V ) = 0}. 

In addition, the attraction domain of the trivial solution is the ball B K , where 
K < H andoo 2 {K) < oo^H). 

Furthermore, the following theorem establishes conditions for global stability: 

Theorem 7 (Theorem 5.8, pp. 79, [33]) Let there exist two functionals 
V(t,<f) and W(t,<f) satisfying the condition of Theorem 6 in any Bh and 

lo\ (|</?(0)|) < V(t, (p), u>i(s) — > oo, s — > oo 
Then the trivial solution of (A.l) is globally uniformly asymptotically stable. 

Using Theorem 6 and Theorem 7 it is straightforward to establish our claim 
of global stability of the trivial solution as stated in Proposition 1. 

rescaled time t — > 2t and rescaled variable y — > y/2 
= -2x(t)+y(t) + 2 7 f[x(t-T/2)] 

(A.2) 

= —Arx(t). 



Proof Consider (3) with 
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Clearly, if the trivial solution of (A. 2) is globally stable, then so is the trivial 
solution of (3). Rescaling was chosen above because it optimizes the bound on 
the stability region. Using (A. 2) and the Liapunov-Krasovskii functional 

V(t, x t , y t ) = + V ^f- + 7 2 f f[x(s)fds, (A.3) 

2 or Jt-r/2 

we obtain 

V(t, z t , y t ) = - {x{tf - 7 2 /kW] 2 } - {lf[x{t - t/2)\ - x(t)f < 0. 

The inequality holds because of the assumptions of Proposition 1, i.e. kf\x(t)\ > 
\f[x{t)}\ (Vrr^O) and | 7 | < kj 1 . 

Condition (1) of Theorem 6 can be satisfied with oui(s) = s 2 /2 and uj 2 (s) = 
[(4r) _1 + t/2]s 2 . Condition (2) is satisfied with cJ 3 (x t ) = V(t,x t ,y t ). We 
note that the set Sips = 0) consists of elements x t = (x t (9) , y t (9)) T G 
C ([-r/2,0],M 2 ) such that f[x(t - r/2)] = x(t) = 0. Define the second func- 
tional as 

W = -x(t)y(t). (A.4) 
Functional (A.4) is bounded in any ball B H and thus condition (3) is satisfied. 

Denote the zeros of / by x* , i.e. f(x*) = 0, and let be a positive constant 
such that, for any of the zeros x* and all x e R, the Lipschitz condition 
\f(x*) — f(x)\ < ki \x* — x\ holds. Define the set E(/i,p) through /i < \\x t \\ < 
H, \x(t)\ < p, and \x* — x(t — r/2)| < p. The last inequality implies that 
\f[x(t — r/2)] | < k L p. In this set we have 

W = 2x(t)y(t) - y{tf - 2 1 f[x(t - r/2)]y(t) + Arx(t) 2 

< 2pH -p 2 + 2^pH + p 2 

< V/2, 

where r < 1/4 and the assumptions of Prop. 1 were used. The last inequality 
follows because p may be chosen such that p 2 + 2(1 + k L /kf)Hp < p 2 /2. Thus, 
W is integrally unbounded as demanded by condition (4). 

Asymptotic stability of the trivial solution of model (A. 2) follows from Theo- 
rem 6. Furthermore, all conditions of Theorem 6 are satisfied for any ball Bh 
and ui(s) — > oo as s — > oo. Thus, the steady state of model (3) is globally 
uniformly asymptotically stable. 
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B Proof of Proposition 2 



Using the identities 



cos(f2r) = -, 



sm(ftr) = sgn(r - — ) : 



cos(2f2r) 



2-b 2 



sin(2f2r) = 2sgn(r — 



and 



4fi 2 - r = -Q sgn(r - ^J)^ 2 ^! + 3^ 
the terms in the denominator of K 2 can be simplified. Define 
D 1 :=D[ + iD'l = b [brtt 2 + (r + Q 2 )e inr] 

then 



D[ = (b 2 n 2 r + r + n 2 ) and D" = sgn(r — —j=)(r + fi 2 )-\/&^--T. 



Applying the identities to 

D 2 := + iD 2 = ^b 2 [2iQb + {4Q 2 - 2iQ - r)e 2inT 
one obtains 

D' 2 = sgn(r - ^)Vb^l(b 2 + 2) - 3fi(6 2 - 2) 



£> 2 ' = 2 (V - 1 + sgn(r - ^)3QVW^j . 
Use of the definitions of D 1 and D 2 allows us to rewrite K 2 (see (39)) as 



+ Im 



6 3 (5 2 ) 2 ft 2 



(D> 1 + iD>>)(D> 2 + iD»)) 

Note that > and therefore the sign of K 2 is given by the expression 
in '{}'-brackets. The definitions of b, b 2 , and 6 3 , i.e. (5) and (10) respectively, 
imply b-b 3 = 7&/'(0)/"'(0) and (6 2 ) 2 = 7£/"(0) 2 . Using this and the definition 



-2b 3 N n 

N n = D[D" + D"D' and C n = T 

D[\D 2 \ 2 



(B.1) 



yields the statement of Proposition 2. 
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C Proof of Proposition 4 



It suffices to determine the sign of N n because the sign of C n is given by the 
sign of the product — bN n (see B.l). For clarity, let us reproduce here (41) 



jV| = (r + Q 2 )(b 4 + 2b 3 + b 2 - 4) + 2nsb 2 (b 3 - 1) 
nix , 



+sgn(r - -^)3fiv / ^ 3 T[(4 - b 2 )(r + fi 2 ) + 2ttsb 2 }. (C.l) 

In writing (C.l), we use the parameterization variable s = Qt > and intro- 
duce the notation TV™ (iV<) for the case r > nn / ' y/r ( r < rmjyfr). Note that 
distinguishing N™ and iV™ allows us to view as function of b, N^=N'^(b, r). 
The main idea behind this proof is to rewrite the sum (C.l) in such a way that 
every term of the sum can be shown to have the same sign. Throughout the 
proof, use is made of the fact that Q and r are strictly positive real numbers 
and that < r < 1/4. We address in turn the four cases of Proposition 4. 

(1) Case r > nn/^/r: 
A) Case b G (1, oo): We need to show that C n < 0, which is true if iV™ > 0. 
Rewrite iV™ as follows: 

AT™ = ( r + Q 2 )(b + 2) [{b - l)(b 2 + b + 2) + 3fi v / ^ 3 T(2 - b) 

+2Qsb 2 {b 3 - 1) + 6Q 2 sb 2 Vb 2 ~^l. (C.2) 

Clearly iV" > for 1 < b < 2. For b > 2 it suffices to show that the 
first term in the sum is positive. However, 

'(6 - l)(b 2 + b + 2) + 3fiv / 6 T3 T(2 - 6)' 



> (b - 2) (b 2 + b + 2 - 30^62^1) (C.3) 
and, using the fact that VL < y^r < 1/2, implies 



(b - 2){b 2 + 6 + 2 - Sfiv 7 ^ 3 !) > (6 - 2)(6 2 - -6 + 2) > 0. 

B) Case 6 G (-oo, -1): To show that iV™ < 0, and therefore C n < 0, we 
rewrite A^> as follows: 

TV" = _( r + Q2)( 6 + 2 )(i _ 6 )( 6 a + 6 + 2 ) _ 2fi 2 r6 2 (l - 6 3 ) 

-Sfiv 7 ^!^ + fi 2 ) (4 - 6 2 ) - QtfrVW^lb 2 . (C.4) 

For — 1 > b > —2 all terms of the sum are less or equal zero and the 
second term is negative definite. To show that iV™ < for b < —2 we 
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make use of the identity 



) )jr _ . 

r + Q 2 = 2Q 2 + sgn(r - —)QVb^l. 



Substituting (C.5) into (C.4) and rearranging we obtain 

- (2s - 12Q)b 2 



N n = b 5 



6 5 + Qtts VW^l 
5 S 4 \b\ 



+ 
+b 3 



-sb 5 + 5Qb 2 + 24Q 2 VV 2 ~^T 

5 



Note that the inequalities 



s > 7T > 3, < -, and — < — < 1 

— — \0\ 

hold for r > nn/y/r with b < —2. Thus, 



(C.5) 



-20fi-4Q 2 \/623i". (c.6) 



(C.7) 



6 5 + 6fisV&^T 9 5 n 

-S — > 71 > 0, 

5 4 \b\ 20 4 

2s-12fi>2(7r-3) >0, 



12. 



17, 



(C.8) 
(C.9) 
(CIO) 



-s6 5 + hVtb 2 + 24fi 2 v / ^ 3 T < — b 5 + —b 2 < 0, 
5 5 2 

where (C.10) holds for b < —2, because 

12 ;5 17 ;2 107 ;5 17, 2n 3 oN n 

—b 5 H b 2 = b 5 H b 2 (b 3 + 8) < 0. 

5 2 80 16 V ; 

Hence, the first three terms of (C.6) are negative definite and this is 
also true for the remaining terms. 

(2) Case r = mxjyfr, n even — > b — 1 : From (C.l) it can be seen immedi- 
ately that C n = N n = 0. 

(3) Case r < mr/y/r, b G (1,2] : We need to show that C n > 0, which is 
equivalent to showing that N™ < 0. Rewrite N™ as 



TV™ = (r + Q 2 ){b + 2) {b - \){b 2 + b + 2)- 3fi(2 - b)^W^l 



— sfl 2 b 2 VV 2 ~^l + 2Qsb 2 



In^/W^i + b 3 - l 



(Cll) 



The expression in brackets in the first term of (Cll) can be bounded from 
above by using the inequalities s > \Jb 2 — 1 < Q < \Jb 2 + 2r — 1, 
r <\, and (r + fi 2 ) < b 2 . It is found that 
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(b - 1)(6 2 + 6 + 2) - 3fi(2 - b)VW^l 

< 4(6 - !)(»» - \ - 1) = 4(6 - 1)0 - ^±i)(6 + ^i). 

Using the above bound, we show that the sum of the first two terms of 
(C.ll) is negative. Consider the first two terms of (C.ll) 



+ Q 2 )(b + 2) \(b - l)(b 2 + b + 2) - 3^(2 - b)VW^l 



-stt 2 b 2 VU^l 
3 



<4(r + Q 2 )(b + 2)(6 - 1)(6 2 - - - 1) - 2tt(6 2 - l^v 7 ^ 7 ! (C.12) 
<0. 



Here, the last inequality follows because both terms in (C.12) are non- 
positive for b G (1, (v / 17+ l)/4] and it holds that 



4(r + tt 2 )(b + 2)(6 - 1)(6 2 - - - 1) - 2tt(6 2 - l)6 2 v / ^T 

Z 



<b 2 (b-l) 



4(6 + 2)(6 2 - - - 1) - 2tt(6 + ljv 7 ^ 7 ! 

Z 



for 6 G ((vl7 + l)/4, 2], where it is straightforward to check that 



4(6 + 2)(6 2 - - - 1) - 2tt(6 + 1)V6 2 - 1 < 0. 

Z 



The remaining term of (C.ll) is also negative on b G (1,2] because 



2Qsb 2 



Ifl^/W^l + b 3 - 1 
o 



< 2Qsb 2 



|(6 2 -l) + 6 3 -l 



2fis6 2 (6- i)(6-2)(6 + -) 



<0. 



(4) Case r < nn/y/r, b G [—2, —1] : By rewriting N™ as 



AT™ = -( r + fi 2 )(6 + 2)(1 - 6)(6 2 + b + 2) - 21M 2 (1 - b 3 ) 

-3nVW^l(A - b 2 ) (r + n 2 ) - 6tt 2 sb 2 Vb^l (C. 13) 



it is immediately obvious that TV" < and therefore C n < 0. 
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